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Abstract 

We show that the Hartle-Hawking vacuum for theories of interacting massive scalars in de Sitter 
space is both perturbatively well-defined and stable in the IR. Correlation functions in this state 
may be computed on the Euclidean section and Wick-rotated to Lorentz-signature. The results are 
manifestly de Sitter-invariant and contain only the familiar UV singularities. More importantly, the 
connected parts of all Lorentz-signature correlators decay at large separations of their arguments. 
Our results apply to all cases in which the free Euclidean vacuum is well defined, including scalars 
with masses belonging to both the complementary and principal series of SO(D, 1). This suggests 
that interacting QFTs in de Sitter - including higher spin fields - are perturbatively IR-stable 
at least when i) the Euclidean vacuum of the zero-coupling theory exists and ii) corresponding 
Lorentz-signature zero-coupling correlators decay at large separations. This work has significant 
overlap with a paper by Stefan Hollands, which is being released simultaneously. 
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I. INTRODUCTION 



While free quantum fields in de Sitter space (dSo) have been well understood for some 
time (see [1] for scalar fields), interacting de Sitter quantum field theory continues to be a 
topic of much discussion. In particular, there has been significant interest in the possibility 
of large infrared (IR) effects in interacting de Sitter quantum field theories [2Tl2"T]. both with 
and without dynamical gravity. 

In [2H] we began to address the specific class of such concerns associated with infra- 
red (IR) divergences of the naive Lorentz-signature de Sitter Feynman diagrams, or more 
generally those concerns that can be addressed in the context of minimally-coupled scalar 
fields with mass M 2 > 0. There we computed one- loop corrections to propagators on 
Euclidean de Sitter (which is just the D-sphere S D ) and analytically continued the results 
to Lorentz-signature. This procedure defines the so-called Hartle-Hawking vacuum of the 
Lorentzian theory [2H] , which on general grounds should be a good quantum state (see section 
|V| . In particular, the analytically continued correlators are expectation values of products 
of operators in a single state as opposed to matrix elements between an "in-vacuum" and a 
potentially different "out- vacuum." We do not attempt to define any notion of S-matrix. 

Because S D is compact, it is a priori clear that Euclidean correlators do not suffer infra-red 
divergences. We showed in [28] that, to one- loop order, the analytically continued Lorentz- 
signature correlators were also finite and decayed at a rate determined by the lightest relevant 
mas a3 The purpose of the current paper is to extend these results to arbitrary iV-point 
functions and to all orders in perturbation theory, again showing that connected correlators 
decay rapidly as the separation between points becomes large. As in [2H], our results will 
apply to all masses for which the free Euclidean de Sitter vacuum is well-defined, i.e. for 
all M 2 > 0, including values in both the complimentary series and the principal series of 
SO(D,l). 

The decay of connected correlators demonstrates that the Hartle-Hawking state is pertur- 
batively stable, and that the Hartle-Hawking vacuum is an attractor state for local operators 
in the sense defined in [28J. To illustrate the main point, let us consider a state constructed 

1 In addition, the one-loop calculations reported in [30] establish that correlators of free-field stress tensors 
decay at large separations. 
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from the Hartle-Hawking vacuum \0)hh with appropriately smeared operators: 

|*>:= f ... / f(Y u ...,Y n )MYi)---MYn)\0)HH. (1) 

Here the Yi are points in dSo, j Y ■ ■ • denotes an integral over de Sitter, and /(Yi, . . . , Y n ) is a 
smearing function which we assume to be supported in a compact domain T>. Now examine 
the correlation function (^>\(f) cr (Xi) ■ ■ ■ o .(Jfjv)| , I r ) with all Xj at large separations from T>. In 
this configuration the correlator is simply a smeared correlation function between 2n opera- 
tors located within D and N operators with large (say, roughly equal) separations \Z\ from 
T> evaluated in the Hartle-Hawking vacuum. Since the associated connected correlators de- 
cay rapidly at large separations, this function approximately factorizes into a product of two 
correlators: one for the points in T> and one for the other points. The former factor is just the 
norm of |tf>, so we have • • • <f><,(X N )\V) ->• (tf |tf > • hh(0\MXi) ' ' ' M^n)\0)hh- 

This means that, as probed by local operators, the excited state becomes indistinguish- 
able from the Hartle-Hawking vacuum. 

We begin by briefly reviewing free de Sitter quantum field theory in section |TT| We then 



address simple tree diagrams in section [HI] which also serves to introduce some useful Mellin- 
Barnes techniques and our choice of (Pauli-Villars) regularization scheme. We address gen- 



eral diagrams in section |IV[ where we establish the desired results for finite Pauli-Villars 
regulator masses (so that all diagrams are finite). Since the infra-red asymptotics are inde- 
pendent of the regulator masses, it is straightforward to take the limit where such regulators 
are removecj^J Some technical material is relegated to the appendices. We close with some 
discussion in section |VJ 

Remark: While paper was being prepared, we received a draft of [33] which reports 
similar results. 



II. FREE DE SITTER QFT 

This brief section serves as a review of scalar quantum field theory in de Sitter and allows 
us to establish our notation. We consider D-dimensional de Sitter space dSn with radius 
2 After subtracting regulator-dependent local counter-terms in order to obtain a finite result. We consider 
theories can be renormalized in this way. One would expect this procedure to be equivalent (up to finite 
local counter-terms) to the renormalization prescription given in |31j . and thus to define a fully covariant 
renormalized quantum field theory in the sense of |32j whenever the flat-space limit is power-counting 
renormalizable. However, we have not analyzedgthis question in detail and save any investigation for 
future work. 



£, which may be defined as the single-sheet hyperboloid in a D + 1-dimensional Minkowski 
space M D+1 . Points on de Sitter satisfy [31] 

VAB X A X B = i\ (2) 

where X A is a vector in the embedding space and t}ab = diag(— 1, 1, . . . , 1) is the usual 
Minkowski metric. Henceforth we will drop the index notation and denote the inner product 
of two embedding space vectors X\ and X2 simply by X\ ■ X 2 . For two points on de Sitter 
located at X\ and X 2 the inner product X\ -X 2 /£ 2 provides a convenient measure of distance 
which we loosely call the embedding distance between X\ and X 2 [I] . The embedding distance 
is related to the length of the chord between X\ and X 2 in the embedding space (with the 
length being proportional to 1 — X\ • X 2 ) and is clearly invariant under the full de Sitter 
isometry group SO(D, 1). The embedding distance satisfies: 

• X 1 ■ X 2 /£ 2 G [— 1, 1) for spacelike separation, 

• Xi ■ X 2 /£ 2 = 1 for null separation, and 

• \X\ ■ X 2 /£ 2 \ > 1 for timelike separation. 

The antipodal point of X\ is simply —Xi\ clearly the embedding distance between antipodal 
points is —1. See Figure 1. 

In this work we restrict attention to massive scalar fields <p a (X). It is convenient to keep 
track of the spacetime dimension with the parameter a = (D — l)/2; the mass parameter a 
is then defined by the equation 

- a(a + 2a) = M 2 £ 2 , (3) 

where M 2 is the bare mass-squared of the field if we assume minimal coupling to the metric. 
There is a redundancy in this definition as ^ is invariant under a — > — (a + 2a); for clarity 
we choose to define a as the positive root 

a-^-a+^-MH 2 ) 112 , (4) 

but all expressions involving a must necessarily be invariant under a — > — (a + 2a). Free 
scalar fields form irreducible representations of the de Sitter group SOq(D, 1) and fall into 
three series 135 1: 
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FIG. 1. The conformal diagram of global de Sitter. The dashed ends are identified. Shown 
are the points X and the corresponding antipodal point —X. Values of the embedding distance 
Z := XY/£ 2 in different regions of de Sitter are labeled; in addition, the dashed red lines denote 
the lightcone with Z = 1 and the dotted green lines denote the lightcone with Z = — 1 [TJ. 

1. complementary series: —a < o < , 

2. principal series: a = —a + ip, p G M., p > , 

3. discrete series: a = 0, 1, 2, . . . . 

We plot <t and — (cr + 2a) as a function of M 2 > in figure 2. Relatively light massive fields 
belong to the complimentary series while heavier fields belong to the principal series. It is 
useful to note that o cc = | — a corresponds to an otherwise massless conformally coupled 
free field. This value lies in the complimentary series so long as D > 2. Discrete series fields 
correspond to massless and tachyonic scalars and we will not consider them here. 

Free massive scalar fields admit a unique de Sitter-invariant Hadamard vacuum |0)f ree , 
commonly referred to as the Euclidean vacuum (it is also the Bunch-Davies vacuum) [U [3] . 
Since the theory is free, the vacuum is completely characterized by its 2-point functions. Let 
us define the function 



a (7 , _ P~ D r (-or) r (a + 2a) 

K[Z) - WY^- r(| + a) 2 1 



o 1 1 + Z 

-a, a + 2a; - + a; 



(5) 



Here 2-^1 ( a > b; c; z) is the Gauss hypergeometric function. In general this function has a 
branch point at Z = 1 and is cut along the positive real axis Z 6 [1, +00). The time-ordered 
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FIG. 2. On-shell values of a and — (a + 2a) in the complex plane for massive scalar fields. The 
red dashed line denotes the path of a for increasing M 2 starting from at a = for M 2 = 0. The 
green doted line shows the path of — [a + 2a) for increasing M 2 starting from — (a + 2a) = —2a 
for M 2 = 0. Relatively light fields with < M 2 £ 2 < a 2 correspond to values of a and — (cr + 2a) 
on the negative real axis and belong to the complementary series. Heavier fields with M 2 £ 2 > a 2 
correspond to complex values of a and — (cr + 2a) on the line defined by Re a = Re (—a — 2a) = —a. 

and Wightman correlation functions of a massive scalar field (f) a {X) are given by 

(OlT^X^pQIO^ee = A^Xi ■ X 2 - ie), (6) 

(0\MXi)MX2)\0)f** = A «(Xi ■ X 2 - ies(X h X 2 )), (7) 

where in <Q the operator ordering is enforced by s(X L ,X 2 ) = +(— ) if Xi is in the future 
(past) of X 2 (see, e.g., [35]). 

At the level of free fields, one may in fact use any member of the one-parameter family 
of 2-point functions found in [H [3] to define a de Sitter-invariant vacuum state. These other 
vacua are usually called Mottola- Allen (MA) or a vacua. However, the non-trivial MA vacua 
do not satisfy the Hadamard or Bunch-Davies criteria; in particular, their 2-point functions 
i) have an additional singularity at antipodal points and ii) have an additional negative 
frequency contribution to the singularity at coincident points [37] . As a result, only the 
Euclidean vacuum extrapolates to the usual Minkowski vacuum in the flat space limit [38J. 
It has been difficult to find consistent extensions of MA vacua to interacting theories - see 
e.g. [36j |39H^3] . For these reasons we will discuss only the Euclidean vacuum in this work. 
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III. SIMPLE TREE DIAGRAMS 



We now proceed to analyze simple connected tree diagrams. As noted in the introduction, 
we compute diagrams on Euclidean S D and analytically continue the results to de Sitter. 
In particular, the Mellin-Barnes techniques used below provide representations of connected 
diagrams Vn{X±, . . . , X n , X N ) on S D in terms of the N(N — l)/2 embedding distances 
Zij = Xi ■ Xj relating the external points. While the are not all independent for general 
N,D, it will often be convenient to use our Mellin-Barnes representation to extend the 
definition of Vn to a function of iV(iV — l)/2 independent variables Z^. The analytic 
continuation can then be performed by analytically continuing in each Zy and evaluating 
= Xi ■ Xj for N points Xj in Lorentz-signature de Sitter space. 

The only subtlety in the analytic continuations will be the presence of branch cuts. As 
noted in section [TTJ for the two-point function this amounts to choosing the appropriate ie 
prescription to construct time-ordered or Wightman correlators, as desired. Much the same 
is true of higher X-point correlators, though the specifics are more complicated to state. 
However, since our only goal is to extract the asymptotics at large Z^, we need not be 
concerned with such details here. The large Z asymptotics are identical on both sides of 
each cut so that all analytic continuations satisfy the fall-off properties derived below. This 
means in particular that our results hold for both Wightman and time-ordered correlators. 



A. The Green's function 

It is convenient for our analysis to use a Mellin-Barnes integral representation of the 
scalar Green's function on S D . Mellin-Barnes representations have proved to be quite useful 
in evaluating Feynman diagrams in flat-space QFT (see, e.g., [44] for an introduction). They 
are especially convenient for deriving asymptotic expansions (see §4.8 of IS]), and it is for 
this reason that we choose to use them here. We review some essential information about 
Mellin-Barnes integrals in Appendix [A] further details can be found an any standard text 
on mathematical methods. 

Starting with the case a < a cc = | — a and a > |, we may write the scalar Green's 
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function 



A a {Z) 



( 47r )a+l/2 r [1 + a + ^ 1 _ a _ a j 



-a + i/, a + 2a + v, —u, 



a — v 



where we use a condensed notation for products and ratios of T-functions: 



r 



01) «2, 

h, b 2 , 



bk 



r (ai) r (a 2 ) • • • r (a^- 
r(6!)r(6 2 )---r(6 fc 



(9) 



or merely r [ai, a 2 , . . . , %] for just a product. In (J8J) the symbol . . . denotes a contour 
integral in the complex v plane. We take as implicit the measure dv/2iri. The contour 
of integration is a straight line parallel to the imaginary axis traversed from — ioo to +ioo 
anywhere within a region called the "fundamental strip" (FS). In general we denote a fun- 
damental strip by its left and right boundaries < l,r >. For the Green's function ([8]) the 
fundamental strip is < a, | — a > which is non-empty due to the restriction a < \ — at. The 
integrand is analytic in v within the FS; beyond the FS it has an infinite number of poles 
due the Gamma functions. By convention we call poles generated by Gamma functions 
T(- • • + v) left poles; likewise, we call poles generated by Gamma functions T(- ■ ■ — v) right 
poles. The fundamental strip is the region between the left and right poles. For this reason 
we do not generally need to write the FS explicitly as it may be inferred from the Gamma 
functions of the integrand. 

The asymptotic behavior of A a (Z) at large \Z\ ^> 1 may be determined by moving the 
contour to the left. The first of the two series of left poles give the leading asymptotic 
terms: 

1 



T {T[-a,a + a}(-2Zy + T[ cr + 2a, — a 



A-n a 

x[l +0 (Z- 2 )]. 



-2Z) 



-<r-2a 



} 



(10) 



The asymptotic behavior for \Z\ near 1 is determined by moving the contour to the right. 
When D is odd, a is an integer greater than or equal to 1 and the leading behavior is given 
by 




l/2-a 



+ r 



\ — a, —a, a + 2a 



a - cr, | + a + a 



x [1 + 0{1-Z)) 



(11) 



When D is even a = | + n, n G No (where No are the non-negative integers) and the two 
sets of poles overlap at v e N yielding double-poles. As a result the pole at v — gives a 
term with logarithmic behavior: 

r(ni fl-Z^~" 



A CT (|^| < 1) 



(47r)n+i 

1 



(4vr; 



ra+1 



[1 + 0{1-Z)\ 

log(l-Z) + 0(l) 



(12) 



1 + a + 2n 
1 + a, 1 + n 

(the first term is omitted when n — 0). 

When cr > <r cc the left-most right pole in ([Sj) lies to the left of the right-most left pole so 
that there are can be no straight contour in between. To arrive at an expression valid for all 
masses, consider again the case a < a cc and move the contour in (pi) to the right past the 



first right pole at v 



a to obtain the expression 



A CT (Z) 



-1 



r 



( 47r )a+l/2 

L 1 

r(«-§) (\ 



-a + v,a + 2a, —u, | 
\ + a + a,\ - a - a 



a — v 



a] 



( 47r )a+i/2 v 2 ; ^ 13 ^ 

In the integral in the first line the contour lies in the interval (max{a, | — a}, min{0, | — a}) 



This interval is non-trivial for a < | — a (since a < 0), and (13) is a valid representation 



of the propagator for any such a . This process can be repeated as needed so that one can 
then increase a as far into the complementary series as desired. The asymptotic properties 



when a > \ — a are again given by ( 10 )-( 12 ). At conformal coupling a 



residue term in (13) survives 



A CC (Z) 



r(a-|) fl-Z 



1/2-a 



a, only the 



(14) 



( 47r )a+l/2 ^ 2 

The behavior of the Green's function at large M 2 3> 1 will be important to our analysis. 
Starting with d8j) we define 



1 



r 



(4tt)«+ 1 /2 

so that the Green's function may be written 

A a (z)= [ ^(f)r(-i/) 



—a + v, a + 2a + v, \— a. — v 
\ + a + a, \ — a — a 



(15) 



1 - Z 



(16) 



At large M 2 ^> 1 the function ip a {v) has the asymptotic behavior 

^2a-\+2v 



I'M 



-r 



a — v 



(1 + 0(M- 2 )) 



(4vr) a + 1 /2 ' 

and as a result the Green's function has the asymptotic behavior 



(17) 



M 



2a-l 



(4tt) q+1 /2 



M 



2v 



(1 + 0(M- 2 )) 



(18) 



Note that (18) contains no left poles; the left poles of the original expression (15) do not 
appear at any finite order in the expansion in inverse powers of M 2 . In the limit M 2 — > oo 
the inequality |M 2 (1 — Z)/2\ > 1 holds for any fixed Z ^ 1, and in this limit the contour in 



(18) may be closed in the left half-plane giving A a (Z ^ 1) = 0(M 4 ). By examining the 



action of (18) integrated against a test function (represented as an MB integral) one may 

6{Z-1) 



determine that (18) is equivalent to 

1 



K(Z) 



1 



M 2 vo\(S 2a ) (1 - Z 2 )"- 1 / 2 
the first few sub-leading terms are 

1 1 5{Z-l) 1 1 



O (M- 4 ) ; 



(19) 



K{Z) 



d 



M 2 vo\(S 2a ) (1 - Z 2 y- X l 2 MUo\(S 2a )dZ 

1 1 d 2 \ 5{Z-1 
+ M®vo\(S 2a )dZ 2 



S(Z-l) 



_(1_Z 2 )«-V2_ 
+ O (M- 8 ) . 



(20) 



(1 - Z 2 )«- 1 /2_ 

Of course, the expansion (20) follows from the fact that the Green's function is the inverse 
of the Klein-Gordon operator using v2 i A/2 = -M" 2 ^^ = -M~ 2 (l + V 2 /M 2 + . . 



B. Pauli-Villars regularization 

Feynman diagrams containing loops in general contain UV divergences which must be 
dealt with through the process of perturbative renormalization. For our purposes it is 
convenient to use Pauli-Villars (PV) renormalization [33]. In PV regularization we replace 
the original scalar Green's function A CT (Z) with the regularized function 

[D/2] 

A^:=M^ + £CiA w (Z). (21) 

8=1 

Here [. . . ] denotes the integer part. This function is nothing more than the original Green's 
functions plus Green's functions of heavy particles with masses M 2 = —pi(pi + 2a). We take 
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the masses Mf to belong to the principal series so that A™ S (Z) will decay for large \Z\ > 1 
at the same rate as A a (Z). The coefficients Cj are bounded functions of the Mf chosen to 
make A^ eg (Z) finite at Z = 1; i.e., to cancel the UV-divergent terms in A a (Z) (including 
the logarithmic divergences that occur for even dimensions). For example, for D = 2, 3 the 
PV-regularized Green's function is 

A^(Z) = A a (Z)-A p (Z), for D = 2, 3 (22) 

while for D = 4, 5 it is 

A^{Z):=A a {Z) + C 1 A pi {Z) + C 2 A Pa (Z), for£> = 4,5 (23) 

where the coefficients satisfy 

d + C 2 = -l, C X M\ + C 2 M\ = —M 2 . (24) 

One may write similar expressions for any dimension (see e.g. [13]) and, if desired, one may 
make further PV subtractions to ensure that A™ g (Z) is differentiable to any desired order at 
Z — 1. Such additional subtractions are useful in dealing with either field-renormalization 
counter-terms or derivatively coupled theories. Below, we assume for simplicity of notation 
that neither of these is present in our theory. However, the analysis is identical in the 
presence of derivative couplings so long as one assumes sufficient PV subtractions to have 
been made to render all diagrams finite at the desired order of perturbation theor}|^J In 
particular, detailed specification of these subtractions is not needed. 

The cancellation of UV singularities has immediate implications for the Mellin-Barnes 
representation of the regulated propagators. Since the short-distance expansion is deter- 
mined by the location of the right-poles, and since right poles with Rev < give terms 
divergent at Z = 1 (where the character of the divergence depends on the location of the 
pole), all such right-poles must cancel; i.e., the fundamental strip for the regularized prop- 
agators may be extended to < a, > without picking up any explicit pole terms of the sort 



that appeared in (13). It follows that for any a < we may write the regularized Green's 
function as 

K cs (z) = [ c g (^)r (-i/) V (25) 



3 For theories that are power-counting renormalizable, one may fix the set of PV subtractions independent 

of the order in perturbation theory. On the other hand, non-renormalizable theories should be treated as 

effective theories. In this case, there is no harm in taking the regularization scheme (i.e., the set of PV 

subtractions) to depend on the order in perturbation theory to which one works. 
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FIG. 3. The single-vertex tree Feynman diagram. 



with 

[D/2] 



i=l 



The function ip™ g (v) is analytic on the interval (Recr, |) in odd dimensions and (Re a, 1) in 
even dimensions. Using the results in appendix [A] one may readily show that the function 
i\)<,{y) - and therefore ip™ g (v) as well - has the asymptotic behavior 

+ iy)\ = e~ 3 ^ 2 \y\- 1+x [1 + Or 1 )] for \y\ » 1. (27) 



Furthermore, the integrand in (25) has only a simple pole at v = which insures that there 
is no logarithmic UV divergence. 

The PV-regularized Green's function A™ g (Z) is a bounded function of Z. Because the 
sphere is compact it follows that using the regularized Green's function to compute corre- 
lation functions yields regularized correlation functions that are bounded functions of the 
embedding distances. The UV divergences of the original perturbation series are recovered 
in the limit Mf — > +oo. We consider theories which can be renormalized by subtracting 
local counter-terms with coefficients depending on the regulator masses M$. As remarked 
in footnote [2] above, one would expect this procedure to be equivalent (up to finite local 
counter-terms) to the renormalization prescription given in |31j, and thus to define a fully 
covariant renormalized quantum field theory in the sense of [32] whenever the flat-space 
limit is power-counting renormalizable. 

C. Single- vertex diagrams 



In this section we compute the connected, single- vertex tree-level Feynman diagram that 
arises lowest order in perturbation theory; see Fig. 3. As stated in section III B for simplicity 
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of notation we assume below that there are derivative couplings. However, the analysis in 
the presence of derivative couplings is essentially identical. 

We find it convenient to first use the PV-regulated Green's functions A^ eg (Z) for our 
computation and then to take the limit where the regulators are removed. While such 
regularization is not in fact necessary for tree diagrams, it has the convenient property 



that it allows us to use the MB representation ( 25 ) which treats all masses uniformly. Our 



discussion below involves a set of fields with mass parameters cij. Note that each a, requires 
its own set of regulator masses M^-, so removing the regulators is the limit M^- — > oo (or 
Pij ->• oo). 

The diagram in Fig. 3 is given by the expression 

V N (X 1} ...,X n ,X N ) = J A£(X! ■¥)■■■ A£(X B • Y)A%{X N ■ Y). (28) 
Here Y is a unit vector and J Y ■ ■ ■ denotes an integral over S D . To compute the right-hand 



side we first expand the Green's functions A^ g (Xj • Y) according to (25): 

1-Xi ■ Y\ Ui 



(29) 



After inserting N copies of this into ( 28 ) the integral over Y becomes 

l-Xx-YY 1 (\ -X n ■ Y\ Vn (\ -X N ■ Y 



This master integral is performed in Appendix [Bj the result is 



(30) 



M 



N 



(4vr) a+1 / 2 f | fl-X 1 -X 2 Y 12 ( \ - x„ ■ :v x 



r [-Pi, . . . , -v n , -vn, 1 + 2a + X) v i] J (a) 



r 



-Oi2, • • • , -a nN , A 1 -u l ,...,A N -u N ,- + a + y^ j u i -y^ / 



a nN 



(31) 



Here L a y-- denotes an integral over N(N — l)/2 integration variables a^-. The a^- are 
labelled according to the corresponding embedding distance Xj • Xj. We use the shorthand 
Ai = XjLi a «r The integration contours lie between their respective left and right poles. 
After performing the shift of variables Vi — Y i/j + A± we obtain 



(Xt) ■ ■ ■ cf> an (X n )<f> aN (X N )) 



('-') 



1 V V \ a 12 / 1 \r V \ a nN 

1 — Ai ■ A.2 \ / 1 — A n • Aat 



r [-a 12 , • • • , -a nN ] V N (a) } (32) 
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with 




V N (a) = (4tt) q+1 / 2 / T [-^, . . . , -v N ] ^{A x +v x ).-. $%{A N + v N ) 



\ + a + J2 Vi + J2 a ij 
l + 2a + y> i + 2Y> ij 



(33) 



Our main task is to determine the fundamental strip of each au variable. The Gamma 



functions in (32) restrict the fundamental strip of each a„ variable to satisfy Reajj < 0. To 



further determine the FS we must determine where the function Vjv(a) ceases to be analytic 
in the a^. When all satisfy Rea^ < the function Vat (a) imposes no further restriction 
on the right side of the fundamental strips. Because of the symmetry of the diagram we 
need only study one variable in detail, say a 12 . As a function of a 12 the function Vjy(a) has 
left poles at 



012 



n, 



012 



o 2 - A' 



n, 



012 



a 



E' 



- n. 



(34) 



In this expression n G N , A[ — Ai — a 12 , etc., and ' a^ = a^ — a\ 2 . We conclude that 
the FS of ai2 is 



a 12 : < max jc^ - A[, a 2 - A' 2 , — i — a — ' a^ j 



> . 



(35) 



Analogous statements hold for the remaining a^. In (35) and below we take the operation 



max to select the greatest real part of any of its arguments. Note in particular that since 
the regulator masses My lie in the principle series (so that Re pij = —a is fixed) the allowed 
strip (35) is independent of the values chosen for the regulator masses My, though it does 



depend on the precise locations chosen for the other contours. 

We can use our knowledge of the fundamental strips of the a^ variables to bound the 
behavior of the diagram at large embedding distances Zij. For example, consider the 
case \Z\ 2 \ ^> 1 and all other Zij ^ 1. We are free to arrange the a^ integration contours such 
that all a^ except a\ 2 are fixed satisfying Rea^ = — e where e is an infinitesimal positive 
constant. In this configuration the FS of a\ 2 becomes 



a 12 : < max {a x , cr 2 } + 0(e), > . 



(36) 



We can therefore move the a± 2 integration contour to a± 2 = max {01,02} + 0(e). In this 
configuration it becomes clear that the diagram decays at least as fast as \Zi 2 \ milx ^ ai ' a ' 2 ^ +0 ^ . 
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More generally we may say that when any embedding distance satisfies \Zij\ ^> 1 the diagram 
decays at least as fast as |Zj J -| <7max+0 W where cr max = max{<7i, . . . , a^} and infinitesimal e > 0. 



The diagram provides the connected part of the PV-regulated N-point correlation 
function (<p ai (Xi) ■ ■ ■ (p an (X n )(f) aN (X N )) to lowest order in perturbation theory Our primary 
goal is to determine the behavior of such connected correlators when the operators are 
taken to large separations, so that several embedding distances become large. From the 
discussion above it follows that the connected PV-regulated correlator decays at least as fast 
as \Z\ arw+0 ^ e \ w here \Z\ is the largest embedding distance between operators. In practice 
the diagram may decay much more rapidly. 



In order to show that the unregulated diagrams have the same IR behavior, we must 
take the limit Mf — y oo where the regulator masses become large. The key step is to recall, 



as noted below (35), that the allowed locations of the contours are independent of the 
regulator masses M^-. We may therefore investigate the large behavior by inserting the 
asymptotic expansion (17) for the ifj pij (Ai + Ui), associated with the propagators for the PV 



regulator masses, into (33) with the a±2 contour fixed at any location allowed by ( |35| ) (and 
analogously for the other ay). To leading order, all dependence on the regulator masses is 
in factors of the form (p 1 'j 2a - 1 + 2A i+ 2u i . The particular power law depends on the location of 
the Vi contours, and the most favorable behavior is obtained by taking the 2^ contours to be 
as far to the left as possible. With this in mind, taking into account certain relevant poles, 
it is straightforward to analyze the large My behavior. The leading term is independent of 
Mij and is obtained by simply replacing every -?// cg with the unregulated ip a ; i.e., just by 
the unregulated expression. Sub-leading terms are suppressed by powers of Mj~ 2 and can 



be neglected. Since the unregulated ip a also satisfy (27) at large imaginary 2/3, the 0(1) 



Mellin-Barnes integral can be analyzed in the usual way to find asymptotic behaviors at 



large \Z^\ dictated by the locations of the contours; i.e., by (35) and its analogues. Thus 
the large \Zij\ behavior of the My — y 00 limit satisfies the same bounds we derived at finite 
M^. In particular, the limiting diagram decays at least as fast as |2'|°"ma X +o(e)^ w h ere \Z\ is 
the largest embedding distance between operators. 
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IV. GENERAL DIAGRAMS 



In this section we analyze connected Feynman diagrams containing loops. We again use 
the PV-regulated propagators of section IIIB| For simplicity of notation we again assume 



that there are no derivative couplings or field-renormalization counter-terms. However, 
the analysis with derivative couplings or field-renormalization counter-terms is essentially 



identical so long as sufficient PV subtractions have been made as described in section III B 



At the technical level, the key step will be to show in section IV B that all diagrams have 
a Mellin-Barnes representation of the following form: 



Vn(X%, • • • , X n , X 



(a) 



r [-a 12 , ■ ■ ■ , -a nN ] V N (a) > , (37) 



N) 

i-x 1 -x 2 y 12 _ _ _ ^ i-x n -x N ^ anN 

where the function Vjv(a) satisfies the following requirements: 

1. Vjv(a) is analytic when all are contained within the region given by the set of 
restrictions 

Re aij e (o- max - Vij(a'), 0]. (38) 

Here cr max is the real part of the mass parameter of the lightest field participating 
in the diagram and Vij(a') is a polynomial function of all the Rea^ variables except 
Re (hence the prime) and has non-negative coefficients. 



2. When the contained in the region (138]) the function Vjv(a) decays at large 

|Imai2| ^ 1 at least as rapidly as 

\V N (x + iy, a 13 , • • • , a nN )\ oc e-^ 2 \y\' 1+x , for \y\ > 1, (39) 

and likewise for the other a^. 

However, let us first discuss the implications of this form and show that it leads to 
exponentially decaying correlators as desired. 

A. Implications of our Mellin-Barnes representation 



To begin, note that the requirement (39) ensures that each integral in (37) converges so 



long as no embedding distance is equal to unity, i.e. when the diagram is evaluated away 

16 



from coincident points. For any = x + iy the integrand in (37) is comparable at large 

\y\ > i to 



-n\y\+iny^y\3/2 



1 — X i ■ Xj 



(40) 



and thus converges absolutely. To evaluate V(X\, . . . , X n , Xn) at coincident points we must 
move some of the contours into the right half-plane. For example, suppose we wish to 
evaluate Vat(Xi, . . . , X n , Xn) at X\ = X 2 . To do so we first move the ai2 contour into the 



right half-plane. In doing so pick up a residue from the pole at a\2 = 0. From (38) it follows 
that V]y(ai2 = 0, . . . ) is regular and so this pole is a simple pole. Upon setting X\ • X 2 — 1 
the remaining contour integral, with ayi (slightly) in the right half-plane vanishes, leaving 
just the residue: 



Vtv(X 2 , X 2 , . • • , X n , X N ) 



(a') 



l-Xo- X, 



1 — Xi ■ Xa 



0-13+^23 



«34 



1 - X, ■ X 



N 



a lN+ a 2N 



X n " X 



N 



a nN 



V N (Q,a 



13, 



a nN , 



r [— ai3, . . . , — a„7v] 
(41) 



Here C » . . . denotes that there is no a\ 2 integral. 



In fact, it turns out that the term on the right-hand side of (41) may be written in form 
(37), i.e. Vn-i(X 2 , ■ ■ ■ , X N ). Said differently, a function Vn+k(Xi, . . . , X N , X N+1 , . . . , X N+K ) 
when evaluated at X N+ \ = ■ ■ ■ = X N+K = Y is itself a function of the form Vn+i(Xi, . . . , X N , Y). 
For example, let us consider when K = 2. Following the procedure outlined above equation 



(41 ) we have 



V N+2 (X U ...,X N ,Y,Y) 



(a') 



1-X 1 -X 2 



X„. ■ X 



N 



a nN 



l- Xl -Y\ ai - N+1+ai - N+2 fl-X N -Y 



a N,N+l+ a N,N+2 



r [—012, • • • , —a-N,N+2] 



V/v+2(oi2, • • • , 0-N,N+2, 0) 



(42) 



In this expression the prime in the (a 1 ) below the integral means that there is no ajv+i,iv+2 



integration. The integrand in (42) is still analytic with respect to the remaining in the 



region given by (38). It follows that after a few cosmetic changes we may write (42) in 
the form of (37). Let us re-label the variables — > c i (here i = 1,...,N), then shift 
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variables ansr+i — > ajjv+i — c$; (42) becomes 



V N+2 (X 1 ,...,X N ,Y,Y) 

| ^ 1-X 1 -X 2 ^ ai2 ^ 1-X n -X N ^ anN ai ' N+1 ^ 1-X N -Y ^ aN ' N+1 

T [-aw, . . . , -a^iv+i] ^(a) J. (43) 
In this expression the integral is over the variables ai 2 , . . . , a^jv+i and V^^(a) is given by 
V^i(a) :— —z 7 / It [c% — ai t N+i, ■ ■ ■ ,cn ~ a N,N+i, —ci, ■ ■ ■ , —cn) 



T [— a^jv+l, • • • , — aAr,7V+l] 7[ c ] 
V/v +2 (ai2, • • • , a n Ar, Ci — Ol,iV+l) • • • 5 C N ~ a N,N+l, Ci, . . . , Cat, 0)}. (44) 

In this expression Jj, . . . denotes contour integration over c\, . . . , cjv- These integrals are 
guaranteed to converge so long as the are within the region for which the integrand of 



(42) is analytic. Although this expression is rather complicated, it is easy to verify that this 
function satisfies requirements (1) and (2) using the asymptotics described in appendix |A| 
The same analysis may be performed for any K > 1 with the same conclusion: the function 



Vn+k(X\, . . . , X N , K, . . . , Y) is of the form of a function Vn+i(X\, . . . , X N , Y) given by (37). 



The last and most important consequence of the form (37) is that the function Vn(Xi, . . . , X n , Xn) 
decays exponentially when evaluated at large embedding distances. For example, suppose 
\Xi ■ X2I 1. A bound on the decay of Vn{X\, . . . ,X n ,X^) can be found in the same 
manner as in the previous section. Let all integration contours except that of a\ 2 be located 



at Rea^- = — e. From (38) it follows that in this configuration 0,12 has a fundamental strip 
at least as large as 

a\i ■ < <W + O(e),0 >, (45) 

so Vn(Xi, . . . , X n , Xn) decays at least as fast as (X\ ■ X2) CTmax+e for any e > 0. 

Furthermore, suppose that removing some vertex results in a disconnected diagram, and 
suppose also that one of the resulting connected components contains none of the original 
external legs. Then this piece contributes only an overall multiplicative constant (which 
is finite at finite regulators masses Mij) to the diagram and does not affect the large Z 
behavior. One may therefore remove such pieces from the diagram when computing cr max 
above. We refer to this process as "trimming," so that the trimmed version of a given 
diagram has all such pieces removed. 
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Obviously, the same result also holds for the other embedding distances. From this 
result it follows that the connected part of a PV-regulated iV-point function - which may be 
described to any order in perturbation theory by diagrams of the form Vn - decays when any 
two operators are taken to be separated by a large distance Z at least as fast as \Z\ ama * + °( e ' , 
where a max is the (real part of the) largest a that appears in any trimmed diagram that 
contributes to the correlator. 



B. Proof of the desired Mellin-Barnes representation 



The proof that all diagrams can be written in the form (37) is through induction. One 



constructs a diagram vertex by vertex, beginning with a single- vertex tree diagram. We have 
already seen that single-vertex diagrams have MB integral representations of the required 
form. Thus one simply needs to show that, upon adding a vertex to an existing diagram 



with the form (37), the new diagram is again of the form (37). We show this below. 



The process of adding a new vertex to an existing diagram is shown schematically in 
Fig 4. Starting with an (N + K)-\egged diagram Vn+k (Xi , . . . , X^, Xjv+i, • • • , Xn + k), one 
attaches a new vertex to the K > 1 external legs Xjv+i, ■ ■ ■ ,Xn + k- One then attaches to 
the new vertex (M — N) new external legs so that the new diagram is an M- legged diagram: 



V M {X U ...,X m ) = J V n+k {Xi> ...,Xn,Y,..., Y)A^ +i (X N+1 ■¥)■■■ K*(X M ■ Y). 



(46) 



This procedure generates all diagrams in which no propagator has both of its ends on 
the same vertex. But adding such one-link loops simply multiplies any diagram by factors 
of A™ g (Y • Y), which are just finite constants due to our PV regularization, and which are 
readily absorbed into the definition of Vm- It thus remains only to show that the diagrams 



generated by the above process satisfy requirements (1) and (2) associated with (37) 



Note that K > 1 in order for the diagram to be connected. Following the discussion in 



section IV A, since K > 1 we know that Vn+k(X\, . . . , X^r, Y, . . . ,Y) can be written in the 



form of some Vat+i(Xi, . . . , Xj^, Y). Inserting this (46) becomes 



V M (X U ...,X_ 



M, 



y 



V N+l (X 1: ...,X N , Y)A%(X N+1 ■ Y) ■ ■ ■ K1A X m ■ Y). 



(47) 
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FIG. 4. The process of adding a new vertex to an existing diagram. 



It is convenient to define n = N — 1 and m = M — l. The integral (46 ) can be computed in 



essentially the same manner as the single- vertex diagram section III We begin by expressing 
both Vjv+i(Xi, . . . , Xn, Y) and the regulated Green's functions in terms of their MB integral 
represent at ions : 



VV+i(Xl, . . . , Xn, Y) 

jT ... jf | ^--Xi-x 2 y i2 ^ i-x n -x N y nN ^ i-x 1 -y j i ^ i-x n -y j n 

r [-ai2, • • • , -a nN , -p±, . . . , -p N ) V N+1 (a 12 , ... , a nN , v x ,...,u N ) \, (48) 



AS*(X, • Y) 



1 - X, ■ Y 



(49) 



Here it is important to keep track of notation. In the first equation we have relabelled 
a^N+i — > v %i i — 1, • • • > N, so that the remaining N(N — l)/2 variables run 012, • • • , a n N- 



In the second expression i runs i = N + 1, . . . , M. Inserting these into (47) we then integrate 



over Y using the master integral Mm (see (31)). After performing a shift of integration 



variables z/j — > Vi + i?j (where = 2^ ? =i fry) we arrive at 



V M (X 1 ,...,X 



(a) J (b) 



1 — X\- Xn + \ 



>1,N+1 



1-Xi-Xs 



1 — X m ■ Xm 



Q12+&12 



1 — X n ■ Xn 



a nN+b n N 



r [— a 12 , . . . , —a n N, —bi2, . . . , —b m M] V^(a, b) 

(50) 
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In this expression the by run over all distinct pairs ij (i.e., 612, . . . , 5 m A/) and 



V M (a,b) 




l + 2a + 5>i + 2D&i. 



Vjv+i(a 12! • • • , a nN , Bi + vx,...,B N + vnW^+i (B n+1 + u N+1 ) ■ ■ ■ ip 1 ^ (B M + v M ) j . 

(51) 



It is now straightforward to determine the region for which the integrand in (50) is 



analytic in the integration variables. The simplest variables to analyse are the by variables 
with N < i, j < M. For these variables the analysis is identical to that performed for the 
single-vertex graph; the result is that the integrand is analytic in the region 



Re b mM e ( max <^ a m - B' a M - B' M , 



a — 



5> 







(52) 



As usual here the prime denotes that b m u is omitted from the sums. For variables with 
1 > i < N and N < j < M one finds 



Re b 1M e (max | cr max - P 1)A r+i(a, b),a M - B' M , -- - a - ^ %j 1 , 



.0 



(53) 



Finally, let us determine the region for which the integrand is analytic with respect to a^- 
while holding the b^ contours with 1 < i,j < N fixed to satisfy Re by = — e. In this 
configuration it is easy to determine that the integrand is analytic when 



Re a X i e (cr max - V 12 (a, b) + 0(e), 0] . 



(54) 



Therefore, we can perform the shift of variables a™ — > Oy — b^ in order to get (50) in the 



form (37). We know that the 6« integrals with 1 < i, j < N will converge in the region given 



by (52)-(54). We see that (52)-(54) satisfy (38), and that all integrals converge sufficiently 



rapidly to satisfy (39). Thus we have shown that Vm(^"i, • • • , Xm) is of the form 37 



C. The regulator limit Mf- — > oo 

Our analysis above is complete at the level of effective theories. In that context, one 
keeps the regulators masses My finite and is careful to ask questions only about physics at 
energy scales much less than Mij. But for renormalizable theories one would like to do more 
and to remove the regulators by sending — > oo before taking the limit of large \Z^\. 
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Such questions are straightforward to address using our Mellin-Barnes representations. 
Note that, as with the tree diagrams discussed in section HI Cj we may study the large 
My limit holding fixed the locations of all contours, subject only to the conditions (38) 



found above. Suppose for the moment that we choose the couplings to be independent of 
the regulators masses My. Then all of the regulator-dependence lies in the functions ip p {v) 
associated with the regulator Green's functions and the coefficients Cj. Note that each term 



in the asymptotic expansion (17) of such functions at large My- again decays exponentially 
away from the real axis (now roug hly as e" 7 ^ 1 ) fast enough for the arguments of sections 
IV A|IV B to holcQ As a result, inserting the expansion (17) into one of our Mellin-Barnes 



integrals (and also expanding the Cj) produces an asymptotic series in the masses My, 
each of whose coefficients is again a Mellin-Barnes integral with the same contours and 
convergence properties as the original expression. 

Of course, the above expansion will in general include positive powers of My- as well 
as negative powers; these are just the expected ultra-violet divergences of the theory. But 
let us suppose that by taking the coupling constants to depend on My in an appropriate 
way the My — > oo limits of correlators become well-defined and finite, at least to some 
fixed order in perturbation theory. This is precisely the assumption that the divergences 
can be cancelled by some set of My-dependent counter-terms. Since coupling constants are 
just overall multiplicative factors in each diagram, it is straightforward to take this extra 
dependence on the My into account. Expanding each coupling in an asymptotic series 
generates a new series, where each term is again a Mellin-Barnes integral of our standard 
form (and with the same placement of the contours). This is true in particular of the term 
that is independent of the My. But this term gives the full My — > oo limit, since all terms 
involving positive powers of My must have cancelled in order to obtain a finite result. The 
usual argument then implies that this term decays as \Z\ ain!lx+0 ^ at large \Z\, where cr max 
is the (real part of the) largest a that appears in any trimmed diagram that contributes to 
the correlator at this order. 



In fact, such arguments require decay only as e ^l/ 2 times an appropriate power law or faster. 
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V. DISCUSSION 



In the above work, we used Mellin-Barnes techniques to determine the asymptotics of 
Pauli-Villars regulated diagrams for massive scalar quantum field theories in de Sitter space. 
We found that connected correlators fall off at large \Z\ at least as fast as does the Green's 
function for the lightest field in the (trimmed) diagram (up to corrections that grow less 
strongly than powers laws; e.g., factors of log \ Z\). Due to the simple way in which changing 
the PV regulator masses interacted with the Mellin-Barnes expressions, it was straightfor- 
ward to show that the same results hold in the Mij — > oo limit in which the regulators 
are removed, independent of the details of any counter-terms required. A similar analy- 
sis using Mellin-Barnes techniques should also be possible in the context of dimensional 
regular ization. 

As described in the introduction, it follows that the interacting Hartle-Hawking vacuum is 
an attractor state in the sense of [25] for local correlators at any order of perturbation theory. 
Our results hold for all masses M 2 > for which a free Euclidean vacuum exists and for 
arbitrary interactions, with non-renormalizable theories being treated as effective theories. 
While for simplicity of notation the calculations were presented only for non-derivative 
couplings, no significant changes are required to analyze derivatively-coupled theories and 
(as usual) derivatives can only strengthen the fall-off at large Z. It would be very interesting 
if our results could be extended to the massless case M 2 = following e.g. the approach of 
|46j . which introduced a new form of perturbation theory on S D . 

Some readers may be concerned by our use of Euclidean techniques. But on general 
grounds the Hartle-Hawking state should be a valid quantum state. In particular, the an- 
alytically continued correlators satisfy the Lorentz-signature Schwinger-Dyson equations. 
Furthermore, the de Sitter analogue jH] of the Osterwalder-Schrader construction implies 
that the Hartle-Hawking state lives in a positive-definite Hilbert space whenever the Eu- 
clidean correlators satisfy reflection-positivity. This in turns holds at least formally whenever 
the Euclidean action is bounded below, and has been rigorously shown in D = 2 dimensions 
for standard kinetic terms and polynomial potentials; see e.g. [IE]. In such cases, it remains 
only to ask how the Hartle-Hawking state relates to other states of interest; e.g, perhaps 
the state defined by the standard in-in perturbation theory in the expanding cosmological 
patch of dSf). This question will be investigated in detail in jl9], where it will be shown 
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that these two states agree for massive scalar fields. 
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Appendix A: Mellin-Barnes integrals 



We write a generic Mellin-Barnes integral as 
f(Z) = 



ai + Aii/, ...,a m + A m u, bi - B u . . . , b n - B n v 
ci + C\v, . . . c p + C v v, d\ - Div, ... ,d q — D q v 



(zy, (Ai) 



where the measure dvjl-ni is implicit and the contour is a straight line parallel to the 
imaginary axis, traversed from —zoo to +zoo, lying between the left and right poles. The 



convergence of the integral (Al ) is governed by the behavior of the integrand at large |Im v\. 



This behavior can be determined from the well-known asymptotic behavior of the Gamma 
function: 



lim r (x + iy) = (27r) 1/2 e-5 |j/| \y\ x ~ 1/2 [l + 0(y- 1 )] . 

[y[->Do 



Let us assume that the all A4, B i: Cj, Di are positive and define 

m n p q 

i=l i=l i=l i=l 

m n p q 



i=l 8=1 i=l i=l 

m n p q 



G = Re 



_i=l i=l i=l i=l 

n p q 

h = u^ A * uw Bi u^~ Ci u^ d 



+ ^(-fn - n + p + q), 



Di 
ij 1 



(A2) 

(A3) 
(A4) 
(A5) 
(A6) 



i=i 



i=l 



i=l 



i=l 



This discussion follows closely the discussion in [50] , 
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and furthermore let Z = Re 1 ® and v = x + iy. With this notation the absolute value of the 
integrand behaves like 

\y\ Fx+G {RHf (A7) 



exp 



-®y- -2~\y\ 



as \y\ — » oo. From this we conclude that the integral (Al) is absolutely convergent when 



1. |$| < En/2. The integral Al defines an analytic function of Z for |argZ| < 
min (it, 4r) . 

2. |$| = Ett/2 and Fx + G < — 1. The integral defines an analytic function for all Z. 



See 1 50 1 for further details. 



Appendix B: Calculation of M.n 



In this appendix we compute the integral 
M N :^M(ux,...,u N ) ' " ' A ' 



l-X n -YY n fl-X N -Y 



"N 



(Bl) 



with n — N — 1. Rather than directly evaluating (Bl) we instead consider the integral 

1-XfY 



A(ai, ...,a n ) :-- 



H YOL r . 



1-X n -Y 



+ 



1-X N -Y 



(B2) 



where «j are arbitrary real parameters and A is a complex number with Re A < 0. The 
quantities A and Ai may be related in a simple way. To do so we use a standard Mellin- 
Barnes formula: 



(At + ■ ■ ■ + A n + A N ) X 



r (-A) J U1 



/■■ 


7 r 


/ 11 1 


Ju„ 



-A + 



, Mi, . . . , 1t n 



i=l 



(^r ■■■(A n )^{A N ) 



(B3) 



Inserting (B3) in (B2) yields 
A(a x , . . . , a n ) 



r (-A) J U1 



"i, 



A + -Ml, • • • , -«n 

i=l 

A< ( «n, A- J. (B4) 
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Written this way Ai is one factor of the Mellin transform of A. 



Let us now return to (B2) and integrate over Y. We use the formula 

•—A roo 



n 



r (-A) J 



d/3/3- 



to write A as 



(2i)~ x f°° 
A{a u ...,a n ) = y c^/T^exp 



^ 1 + ai 



,+i/3VY 



(B5) 



(B6) 



Y 



where V = a±X\ + • • • + a n X n + Xjy. The integral over Y can be written in terms of the 
Bessel function: 



-i/3VY 



(27T) 



The Bessel function may be written as a Mellin-Barnes integral 



r 



2, 



(B7) 



(B8) 



inserting (B8) into (B7) yields 



V J fi 



1 + a + fi 



2t/2\ M 



4 



(B9) 



After inserting (|B9|) into (]B6|) we may integrate over using the inverse of (B5) 

(77 - iO) A . 



d/3(3' 



(BIO) 



Convergence of the integral over (3 requires Re (A — 2/i) < 0. The result is 



A{a ly ...,a r 



21-A^q+I 

TFaT 



—fx, 2/i — A 
1 + a + fi 



\-2fi 



i=l 



Next we perform a number of manipulations in order to tidy up (Bll). First note that 



V 2 = a\ H h a\ + 1 + aia 2 AiX 2 H h a n I n Xjv 

n \ 2 

1 + ai ) + Zai^C^Xa - 1) + ■ ■ ■ + 2a n (X n X w - 1). 
i=i / 

It is convenient to use IB3I to write 

[2« 1 a 2 (X 1 -X 2 -!) + ■■■ + 2a n {X n X N - 1)] 



(B12) 



a,- 



(B13) 
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Inserting this into (Bll) yields 
A(ai, ...,ot n 



21-a^q+i 



r(-A) J, 




r 



2/i — A, —[/, + w, —w 
l + at + fi 



\-2w 



(2i 



-2/1 



a: 



1=1 



[2cna2(XiX a - 1) + • • • + 2« n (X n X iV - l)] w j (B14) 

We can now integrate over \i. First we use the Gamma function duplication formula 

x, x + - = 2 1 ~ 2x ^T (2s) (B15) 

on the Gamma function T (2/i — A); second we use the Gauss summation formula [50J written 
here as a Mellin-Barnes integral: 



a + p,, b + p,, d — fi 

C + fJL 



a + d, b + d, c — a — b — d 
c — a, c — b 



(B16) 



valid for Re (c — a — b — d) > 0. Cleaning up we have 



A(ai 



j ... : QL n 



21+20^0+1/2 

r[-A, l + 2a + A] J ,„ 

A-2to 



r 



2w — A, - + a + A — w, —w 



i=l 



1 — X x • X 2 \ / 1 — X n X N 

«i«2 I H ha n 



(B17) 



The next series of steps is simple but rather cumbersome to transcribe. We expand both 



the term in parentheses and the term in square brackets in (B17) using the Mellin-Barnes 



expansion (B3 ). Within the parentheses there are n+1 terms, so the Mellin-Barnes expansion 



of this quantity has n integrations. Likewise, the term in square brackets has N(N — l)/2 
terms so the Mellin-Barnes expansion of this quantity has N(N — 3)/2 integrations. After 
performing some shifts in the integration variables (taking care not to shift a contour through 
a pole) and relabelling we obtain the following expression: 

A{ai, ...,«„) = 

21+20^0+1/2 



r[-A, i + 2« + a] y w 



It*! 



1 - X x ■ X 7 



1 — X n ■ Xn 



T[-h 



12, • • • , —h n N 



hi - h ni - /jL n , h m - a + \ + a + a - || 



(B18) 
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In this expression there is a total of n integration variables . . . ,/i n and N(N — l)/2 
variables hy. The latter are labelled such that each factor of (1 — Xi ■ Xj)/2 is raised to the 
power hij. 

The convergence of each Mellin-Barnes integral may be evaluated using the technique 
described in Appendix [Al Each integral converges absolutely for all (1 — Xi ■ Xj)/2 ^ 1. 



The expression (B18) defines a single- valued function of the inner products Xi ■ Xj for all 



complex values of Xi ■ Xj away from the cuts Xi ■ Xj £ [1, oo) 



Both (B4) and (B18) equate A with an n-fold Mellin transform with parameters «j. It is 



easy to see that the integration contours of the two expressions - those of the Ui in the former 
expression and \ii in the latter expression - may be taken to be traversed in the same places 
in their respective complex planes. Now recall that the Mellin inversion theorem states that 
for a given choice of integration contour the Mellin transform of a function is unique [50J. 
It follows that we may identify the integrands and equate U\ = . . . , u n = // n . The final 
step is to relabel 



Hi 



Ui -> Ui, for i = 1, . . . , n, A -> v N + ^ v h a^, 



(B19) 



i=i 



which yields 

M(v 1 ,...,v N ) 



(4?r 



,a+l/2 



T -X 1 -X 2 X ai2 



r [-ui, . . . , -u N , 1 + 2a + J2 J[a) 

r 



1 — X n ■ Xn 



a nN 



-a 12 , . . . , -a nN , Ax - ui, . . . , A N - v N , - + a + ^ Vj - ^ 
In this expression Ai := Y^j=i a ij- 



(B20) 
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